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1.  INTRODUCTION 


This  report  describes  a  project  that  investigated  the  field  distributions  produced  by 
incident  electromagnetic  waves  near  dielectric  structures.  The  structures  are  flat  slabs, 
hollow  slabs,  hollow  cylinders,  and  hollow  cones,  all  with  finite  dimensions.  The 
structures  are  idealized  but  are  realistic  models  for  analyzing  the  electromagnetic 
effects  that  occur  near  practical  radome  configurations.  Because  wave  polarization 
direction  strongly  influences  radome  boresight  error,  the  project  emphasizes  the 
polarization  dependence  of  wavefront  distortions.  The  work  in  the  project  dealt  with 
two  approximations  commonly  made  in  analyzing  radomes:  (1)  treat  the  radome  as 
locally  flat  by  approximating  its  transmittance  as  that  of  an  assembly  of  infinitely  broad 
dielectric  sheets;  (2)  omit  guided  waves. 

Radomes  produce  several  electromagnetic  effects.  The  simplest  are  reflection  and 
refraction,  which  involve  essentially  plane  waves.  In  addition,  radome  edges  or  vertices 
scatter  non-plane  waves.  Finally,  radomes  generate  guided  waves,  with  propagation 
constants  that  differ  from  those  of  free-space  waves.  These  effects  depend  on  the 
radome  configuration,  and  the  frequency,  polarization,  and  direction  of  the  incident 
wave.  For  a  coherent  externally  incident  wave,  the  distinct  radome-generated  waves  are 
coherent;  they  interfere  to  produce  a  complicated  field  distribution  near  the  radome. 
The  practical  consequences  are  boresight  error,  attenuation,  and  increases  in  sidelobe 
levels.  Consequently,  radomes  strongly  influence  missile  guidance.  Moreover,  radome 
effects  depend  on  temperature,  which  increases  nonuniformly  in  flight. 

The  effects  of  radomes  on  electromagnetic  waves  motivated  this  research.  The  project 
had  two  goals.  One  was  to  provide  a  physical  basis  for  numerical  radome  analysis  that 
describes  the  complicated  fields  near  radomes  in  terms  of  simpler  constituent  waves, 
analogous  to  geometric  diffraction  theory,  which  has  been  developed  mainly  for  metallic 
scatterers.  The  other  goal  was  to  develop  a  new  guidance  concept  that  utilizes  antennas 
embedded  in  dielectrics;  this  approach  is  analogous  to  integrated  optics.  The  two  goals 
are  closely  related  because  one  seeks  to  understand  guided  waves  on  finite  structures  and 
the  other  seeks  to  exploit  the  guided  waves. 

Our  approach  involved  theory,  computation,  and  experiment.  The  theory  includes  a 
diffraction  theory,  with  integration  over  a  Huygens  surface,  and  a  generalization  of  the 
moment  method  which  has  been  developed  by  J.  Richmond  for  dielectrics.  The 
experiments  involved  probing  the  complex-valued  field  in  a  radome-bounded  region  with 
a  small  antenna  to  evaluate  transmittance  and  identify  constituent  fields  from  their 
interference  patterns,  and  then  comparing  measured  and  computed  values. 

Section  2  of  this  report  describes  radome  boresight  error  and  surveys  some  analytical 
methods  that  are  approximate.  It  describes  the  validity  of  the  methods  and  describes 
omissions  in  them;  specifically,  guided  waves  and  a  local  flat  sheet  approximation  are 
significant.  Polarization  also  is  a  variable.  This  section  describes  an  integrated 
antenna-radome  concept  that  exploits  guided  waves  and  refracted  waves  on  dielectric 
slabs. 


Section  2  also  gives  a  new  moment  method  theory  for  hollow  curved  shells,  like  cones. 
This  theory  is  based  on  Richmond's  vector  and  potental  formulation,  but  it  is  new 
because  it  develops  closed  form  expression  for  matrix  elements  of  doubly-curved  hollow 
shells  rather  than  cylinders. 

To  connect  moment  method  theory  and  radome  problems,  Section  3  first  gives  numerical 
results  for  slabs  and  hollow  wedges.  Analysis  of  slab  data  shows  that  guided  waves  are 
excited.  Analysis  of  wedge  data  predicts  boresight  error  properties,  especially  polariza¬ 
tion  dependence  of  the  algebraic  sign  of  error.  This  polarization  dependence  was 
observed  in  measurements  made  for  an  axially  symmetric,  pointed  radome;  data  are 
given  in  paragraph  3.4.  In  addition  for  plane  wave  incidence,  wavefront  phase  and 
intensity  were  measured  in  the  region  bounded  by  the  axially  symmetric  shell.  The 
measured  wavefront  data  were  compared  with  values  computed  by  a  surface  integration 
method.  Discrepancies  between  these  calculations  and  measurements  were  shown  to  be 
caused  by  omitting  a  guided  wave. 

Section  4  describes  measurements  and  calculations  for  fields  behind  a  hollow  finite 
cylinder.  The  new  moment  method  was  used.  Two  formulations  were  compared.  One 
included  only  radial  field  components.  The  other  included  axial  and  radial.  Dis¬ 
crepancies  between  calculation  and  measurements  were  smaller  when  both  radial  and 
axial  components  were  used. 

Section  5  described  experiments  and  moment  method  results  for  a  hollow  cone. 
Accuracy  was  greater  when  both  axial  and  radial  fields  were  used. 
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2.  ANALYSIS 


Radomes  reduce  missile  guidance  accuracy  by  causing  boresight  error,  which  is  an  error 
in  the  direction  of  a  wave  source,  either  an  emitter  or  a  reflector.  Although  progress 
has  been  made  in  analyzing  boresight  error,  most  methods  are  approximate,  so  they  are 
accurate  for  some  configurations  but  not  for  others.  Unfortunately  extensive  tests  of 
accuracy  are  lacking. 

Relatively  simple  methods  are  adequate  if  the  radome-endosed  antenna  is  an  order  of 
magnitude  larger  than  the  wavelength  (References  1  and  2).  For  example,  consider  the 
direct  ray  method  shown  in  Figure  2-1.  Propagation  to  an  observation  point  in  the 
radome-bounded  region  is  described  by  a  single  ray  through  the  point  and  the  radome 
where  the  ray  intersects  the  radome.  Transmittance  for  the  ray  is  that  of  a  flat  sheet 
that  depends  on  position  and  incidence  direction.  The  transmittance  is  the  ratio  of  the 
complex-valued  field  amplitude  with  a  dielectric  structure  present  to  that  of  the 
incident  field. 

For  antennas  with  a  diameter  of  about  five  wavelengths  the  direct  ray  method  is 
inadequate.  Accurate  boresight  error  was  calculated,  however,  by  integrating  over  an 
incident  wavefront,  a  Huygens  surface.  Again  the  transmittance  was  that  of  flat  sheets, 
but  for  converging  rays,  as  in  Figure  2-2.  For  antennas  with  apertures  approximately  a 
wavelength  in  size,  like  broadband  antennas,  surface  integration  gives  marginal 
accuracy.  See  Reference  3. 


AN  0060  A 


Figure  2-1.  The  Direct  Ray  Method.  A  single  incident  wave  normal 
(or  ray)  is  associated  with  a  point  on  the  receiving  aperture. 

The  radome  is  a  surface  of  revolution,  but  or  y 
a  plane  is  sketchec 


Figure  2-2.  Surface  Integration  Method.  The  radome  is  a  surface  of  revolution. 

Integration  is  over  a  portion  of  the  incident  wavefront. 

To  evaluate  the  accuracy  of  analytical  methods,  measured  and  computed  boresight  errors 
can  be  compared,  but  the  enclosed  antenna  integrates  over  its  extent  and  obscures 
transmittance.  Moreover,  experiments  with  specific  antennas  lack  generality.  The 
influence  of  aperture  size  can  be  minimized  however,  by  probing  the  wavefront  in  the 
radome-bounded  region  with  a  small  antenna.  Proving  gives  explicit  data  on  the 
transmittance  over  a  small  region  and  directly  tests  the  propagation  description  and  its 
approximations.  See  References  2  and  3. 

The  transmittance  depends  on  the  wave  mechanisms  present.  Direct  rays  and  converging 
rays  correspond  to  plane  waves  with  free  space  propagation  constants.  Edges  and 
vertices  give  approximately  cylindrical  and  spherical  waves,  again  with  free  space 
propagation  constants.  In  addition,  dielectric  structures  support  guided  waves.  The 
influence  of  guided  wave  amplitudes,  which  decrease  exponentially  with  distance  from  an 
infinite  slab,  are  slight  for  antennas  that  are  reasonably  large.  However,  guided  waves 
seem  significant  for  small,  broadband  antennas.  See  Reference  3.  In  addition,  guided 
waves  have  potential  for  new  systems  with  antennas  embedded  in  radomes,  in  analogy  to 
integrated  optics. 
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2-1  GUIDED  WAVES;  A  NEW  APPROACH  TO  GUIDANCE  AND  DIRECTION  FINDING 

The  new  approach  consists  of  an  array  of  self-referencing  interferometers  (References  3 
and  4).  A  self-referencing  interferometer  consists  of  a  dielectric  slab  and  a  set  of 
antenna  elements  distributed  over  the  slab.  The  interferometer  operates  because  a  wave 
incident  on  a  slab  excites  guided  waves  as  well  as  reflected  and  refracted  waves,  as 
shown  in  Figure  2.1-1.  The  guided  waves  have  propagation  constants  that  differ  from 
those  of  the  other  waves,  so  an  interference  fringe  pattern  is  formed  (Figure  2.1-21.  The 
maxima  of  intensity  are  chosen  as  locations  for  sensing  antennas  because  the  gain  of  the 
antenna  is  an  order  of  magnitude  greater  than  that  of  the  antenna  in  free  space. 
Because  the  location  of  maxima  depends  on  frequency,  several  distinct  antennas,  each 
for  limited  frequency  ranges,  can  be  utilized. 
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5S  =  Ztt/IKj  -  K„  sin  a) 


Figure  2.1-2.  Intensity  Fringes  Are  Formed  Near  a  Slab 
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Figure  2.1-3.  Antennas  Embedded  in  a  Nonuniformly  Thick  Slab 

In  an  earlier  project,  we  analyzed  the  fringe  patterns  for  flat  slabs  with  the  moment 
method.  We  fabricated  an  interferometer  on  a  slab  (Figure  2.1-3)  and  measured  the 
farfield  radiation  pattern  (Figure  2.1-4).  The  computed  pattern  is  shown  in  Figure  2.1-5. 
We  also  combined  two  interferometers  to  form  a  hollow  wedge  and  measured  the  sum 
and  difference  patterns  for  a  monopulse  arrangement  (Figures  2.1-6  and  2.1-7).  The 
deep  minimum  in  the  difference  pattern  is  in  the  plane  of  the  slab.  For  an  axially 
symmetric  structure,  interferometers  would  be  arranged  as  in  Figure  2.1-8. 

The  earlier  project  emphasized  transverse  electric  fields  or  perpendicular  polarization. 
In  this  project  we  extended  analysis  to  include  parallel  polarization  for  transverse 
magnetic  waves  on  slabs  and  investigate  axially  symmetric  radomes. 

2.2  MOMENT  METHOD 

The  moment  method  procedure  for  slabs  and  hollow  wedges  was  described  in  General 
Dynamics  Electronics  Division  Report  R-81-125  (Reference  5);  it  is  based  on  a  theory 
originated  by  Richmond  (References  6  and  7)  wherein  a  finite  slab  is  represented  by  a 
finite  number  of  infinitely  long  cylinders  called  cells.  An  integral  equation  represents 
the  total  field  E^,  where  E^  is  the  sum  of  the  incident  and  scattered  field  E1  and  E^, 


ET  =  E1  +  Es 


(2-1) 


The  scattered  field,  from  a  vector  potential  A  and  a  scalar  potential#,  is 

Es  =  i<wA  -V'#,  (2-2) 
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Figure  2.1-4.  Intensity  for  Two,  Colinear  Dipoles  in  Nonuniform  Slab. 

Frequency  9.4  GHz. 
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Figure  2.1-5.  Computed  Field  in  a  Quarter-Inch  Tapered  Slab  with  Wires  and 

Dipole,  at  Location  of  Dipole 


AMR-015 

OFF  AXIS  ANGLE  (DEG) 


Figure  2.1-6.  Difference  Mode  Pattern  for  Two  Interferometers  Formed  into  a  Wedge. 
Frequency  9.4  GHz.  The  off-axis  angle  is  measured  from  the  symmetry  axis  of  the  wedge. 


Figure  2.1-7.  Sum  Mode  Pattern  Formed  by  Two  Interferometers  Formed  into  a  Wedge. 
Each  slab  tapered  thickness  and  had  two  colinear  dipoles  on  its  surface. 
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(a)  INTERFEROMETERS  ON  MISSILE  NOSE 


(b)  TWO  INTERFEROMETERS  (cl  THREE  INTERFEROMETERS  (d)  FOUR  INTERFEROMETERS 


AMR  01  7 

Figure  2.1-8.  Installations  of  Self-Referencing  Interferometers  on 
Missile  Nose.  The  interferometers  in  (b)  can  be  diametrically 
opposite  for  a  rolling  airframe. 

where  i  is  J~-\,  <o  is  the  radian  frequency,  and  time  dependence  is  exp(-icat).  The  prime 
indicates  that  the  gradient  is  taken  at  the  observation  point.  The  vector  potential  is 

4  ;rA  =  -i(*i0E0  / («-l)g  Et  dv,  (2-3) 

where  k  is  dielectric  constant,  g  is  r_1  exp(ikr),  k  equals  2tt/X,  A  is  wavelength,  and  r  is 
the  distance  between  the  integration  and  observation  points.  The  scalar  potential  is 

4  7T$  =  -  J* gV.(K-l)  Et  dv  (2-4) 

Equation  (2-1)  is  an  integral  equation.  It  is  reduced  to  a  set  of  simultaneous  equations  by 
decomposing  the  scatterer  into  elementary  cells  that  are  small  enough  to  justify 
assuming  eT  constant  in  each  cell.  From  Equations  (2-1)  and  (2-2)  we  obtain  the  set  of 
simultaneous  equations,  as  follows.  For  a  typical  cell,  which  is  labeled  with  index  m,  we 
enforce  the  integral  equation.  Therefore,  at  the  center  of  each  cell. 


(2-5) 
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The  index  m  ranges  from  1  through  N,  where  N  is  the  number  of  cells.  Em  is  a  sum 
because  all  cells  contribute  to  the  field  at  any  one  cell.  Therefore,  we  revise  the 
notation,  so 


ET  -  ES  =  E1  (2-6) 

m  mn  m 

The  integrals  for  Equation  (2-2)  were  evaluated  in  Reference  5  on  the  assumption  that 
E^  is  transverse.  In  the  coordinate  system  of  Figure  2.2-1,  the  assumption  for 
transverse  fields  is  that  the  longitudinal  or  z  component  Ez  is  zero.  This  section 
generalizes  the  analysis  of  Reference  5,  giving  the  fields  that  result  when  both  Ex  and  Ez 
are  included  and  assumed  nonzero.  We  present  only  the  results  here  since  the  basic 
derivations  are  contained  in  Reference  5.  The  derivations  with  Ez  not  zero  are  similar 
but  somewhat  more  involved,  and  the  gradient  operation  introduces  additional  terms, 
which  can  be  described  as  a  coupling  of  field  components. 

We  consider  dielectric  rings  as  in  Figure  2.2-1.  A  set  of  rings  with  distinct  radii 
generates  a  hollow  cone,  and  a  set  of  rings  with  equal  radii  generates  a  finite  cylinder. 
The  coordinate  system  and  terminology  of  Figure  2.2-2  are  used  in  the  following 
discussion. 
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Figure  2.2-1.  Dielectric  Ring  and  Coordinate  System 
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Figure  2.2-2.  Coordinates  for  a  Cell.  A  typical  Cell  has  Dimensions: 

Radial  2 8p,  Circumferential  25<i>,  and  Longitudinal  2fiz 

2.2.1  SCATTERED  FIELD  FOR  DIAGONAL  ELEMENTS 

2.2.1. 1  VECTOR  POTENTIAL  -  The  scattered  field  at  cell  m  produced  by  the  vector 
potential  within  cell  m  is 

EAm  =  2(«-l)  ji  Jk5z5ppm«<i> - 8p8<t>s* 

Am  (  12  7T  8  z2  2 

+  ~  -  log  (S2/2  8  z)2  j  |  *2"7* 

where 

s2  =  8f>2  + 

2 


Equation  (2-7)  is  Equation  (2-39)  of  Reference  5.  The  expression  for  the  diagonal  terms 
is  identical  to  that  for  Ez  =  0.  Of  course,  Equation  (2-7)  is  a  vector  equation.  That  is. 
Equation  (2-7)  applies  to  Ex  and  Ejr. 
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2.2. 1.2  SCALAR  POTENTIAL 


-  The  x  and  y  components  of  the  scattered  field  are 


xm 


ym 


( 1 — #c ) 

I  2 


1  +  ( 1  -In  q|)  <1  -q2)  +  q3 


cos 


fi 


m 


I  2 

\2pmj 

(1  -In  qj)  q3 


sin2</>m  '  Et 


|  'xm 


4i 


-  ( « - 1 )  ts  8pS</>cos</>m  et  , 


zm 


4i 


-  ( k  - 1 )  ts  8p8<(>  sin  «f»m  E 


zm 


(2-8) 

(2-9) 


where  ts  is  8p/coscv,  ais  the  half  angle  of  the  cone  if  the  scatterer  is  a  cone,  qj  is 

(s/2  8z)2,  q2  is  |  (s/pm  8<i>)2,  and  q3  is  (8ppm  8<j>)/3n8  z2.  Equation  (2-8)  states  that 

Es  is  affected  by  ET  . 
xm  zm 

2.2.2  SCATTERED  FIELD  FOR  OFF-DIAGONAL  ELEMENTS 

2.2.2. 1  VECTOR  POTENTIAL  -  The  vector  potential  of  cell  n  generates  a  scattered 
field  x  component  at  cell  m  as  follows: 


'Axm 


=  ^  (k-1)  k28p8<(>  8  z  pn  -  sine  (kb8pB-l)  sine 


ikB 


kB  ^  (z^~Zfyj)8z 


E  (2-10) 
xn 


where  sine  z  is  (z-^  sinz),  b  is 


f’n"(Jm  cos  ^n-</W 


,  and 


B  = 


P  n  +  2pnpm  cos  (zn-zm)  +  (zn-zm) 


1/2 


(2-11) 


T  T 

The  z  component,  due  to  A,  is  also  given  by  Equation  (2-10)  with  E  replaced  by  E  . 

xn  zn 

2. 2. 2. 2  SCALAR  POTENTIAL  -  The  scalar  potential  generates  an  x  component  at 

cell  m  due  to  the  nth  cell  as  follows: 


Es  =  =  (K-l)8p8<i>8  z  cos<in  ei,<B  B-2x 

xm  17  n 


2-10 


2(B'l-ik)c  +  ikpr 


(2B  1-ik)B'1  be  cos^ 


1  ET 

{  xn 


2i  I 

+  («-l)  -jp-  Sp84>  ts  kpn  eikB  B'2  zmn(ik-2B-1)B~1  e  cos</>m 


I 


t-  zmn  (ik-2B_1)PnB_1  sin  (c6n-{i>m)  sin<i>n  j  E 

where  c  is  p;ncos</>m  -pncos</>n,  and  zmn  is  zm-zn. 

The  z  component  is 


zn 


ES  =  (K-D^SpS^tgkpn  eikB  B-2 


zm 


1  —  (ik—  2 B  l)zmn  znm  B  * 


jzn 


+  (k-1)  z  cose/>n  j-ikpnb  (2B  l-ik)B"l-(B-1-ik)  |  znm  e'kBB~2  E 


(2-12) 


(2-13) 
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3.  DIELECTRIC  SLAB,  WEDGE,  AND  SHELL 


The  wedge  and  slabs  described  in  the  following  sections  were  analyzed  because  they  have 
no  curvature  and  introduce  effects  caused  by  finite  structures.  These  effects  are  guided 
waves,  edge-scattered  waves,  and  polarization  dependence.  Wedges  can  be  analyzed  by  a 
moment  method  that  is  rigorous  and  has  been  extensively  tested  by  comparison  with 
experiment.  Because  the  flat  slabs  and  flat  wedge  wells  lack  circumferential  curvature, 
polarization  effects  are  separate  for  transverse  electric  (TE)  and  transverse 
magnetic  (TM)  waves. 

In  1980  and  1981,  we  analyzed  fields  near  a  finite  slab,  length  10  inch  and  thickness 
1/8  inch  for  wavelength  1.259  inch  (Reference  5).  The  total  and  scattered  fields  were 
computed  by  the  moment  method  for  grazing  incidence  but  only  for  perpendicular 
polarization.  The  main  result  was  a  resolution  of  the  field  distribution  into  constituent 
fields,  which  were  the  incident  wave,  edge-scattered  cylindrical  waves,  a  guided  wave 
propagating  in  the  direction  of  the  incident  wave,  and  a  guided  wave  propagating 
backward  toward  the  source. 

However,  practical  radomes  usually  resemble  axially  symmetric  shells  more  closely  then 
wedges.  Therefore,  we  also  investigated  the  fields  near  an  axially  symmetric  shell.  The 
purpose  in  investigating  an  axially  symmetric  shell  was  to  determine  whether  the  fields 
can  be  decomposed  into  simple,  constituent  waves  that  resemble  those  on  slabs  and 
wedges.  This  decomposition  is  analogous  to  solving  canonical  problems  in  the  geometric 
theory  of  diffraction. 

3.1  HOMOGENEOUS  DIELECTRIC  SLAB,  LENGTH  3.18  WAVELENGTHS 

The  moment  method  was  utilized  to  compute  the  fields  near  a  slab  with  the  parameters 
given  in  Table  3.1-1.  Slab  coordinates  are  defined  in  Figure  3.1-1.  Incidence  was 
grazing,  and  thr.  wavelength  used  was  1.259  inch.  Calculations  were  done  for  both 
parallel  and  perpendicular  polarizations,  which  correspond  respectively  to  transverse 
magnetic  (TM)  and  transverse  electric  (TE)  fields.  The  slab  was  divided  into  33  cells. 


Table  3.1-1.  Parameters  of  Homogeneous  Slab 


Length: 

4  in. 

Thickness: 

0.125  in. 

Dielectric  constant: 

2.6  in. 

The  magnitude  of  the  total  field  at  the  cell  centers  is  shown  in  Figure  3.1-2  for  parallel 
and  perpendicular  polarizations.  ET  is  larger  for  perpendicular  polarization.  The  rapid 
oscillations  for  perpendicular  polarization  can  be  explained  by  interference  between 
guided  waves  propagating  in  opposite  directions. 
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Figure  3.1-1.  Coordinates  for  a  Slab 

The  phase  values  of  are  shown  in  Figure  3.1-3.  The  phase  increment  with  position  is 
obtained  from  the  propagation  constants  kg  of  waves  guided  by  an  infinitely  broad  slab. 
Transcendental  equations  are  solved  which  have  the  following  form  for  even  TE  modes 
(Reference  8): 


-  (ua/2)  cotan(ua/2)  =  (va/2) 


(3-1) 


u2  +  v2 


k2(K-  1) 
o 


(3-2) 


where  a  is  slab  thickness,  kQ  is  2rr/A,  A  is  the  wavelength,  and  k  is  dielectric  constant. 
For  odd  TM  modes  Equation  (3-1)  is  replaced  by 

(ua/2)  tan(ua/2)  =  (va/2)  (3-3) 


The  propagation  constants  for  the  slab  of  Table  3.1-1  are  given  in  Table  3.1-2. 

Table  3.1-2.  Guided  Wave  Propagation  Constants  for  Slab  of  Table  3.1-1 


Mode 

kz 

Lowest  even  TE 

1.758  7r/inch 

Lowest  even  TM 

1.634  7r/inch 

For  the  TE  mode,  kz  gives  phase  delay  316*/inch,  compared  to  the  moment  method 
results  of  308*/inch  (Figure  3.1-3).  For  the  TM  mode,  kz  gives  294Vinch,  compared  to 
the  moment  method  which  gives  290*/inch.  This  close  agreement  supports  the  inter¬ 
pretation  of  slab-guided  modes  even  in  the  short  slab  of  length  3.18  wavelength.  The 
difference  between  the  phase  delay  for  perpendicular  and  parallel  polarizations  will  be 
shown  to  be  significant  for  axially  symmetric  radomes. 
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Figure  3.1-2.  Total  Field  Magnitude  Inside  4-Inch  Long,  0.125-Inch  Thick  Slab; 
Dielectric  Constant  2.6;  Frequency  9.40  GHz;  Grazing  Incidence 


The  distribution  of  the  total  field  near  the  slab  is  shown  in  Figure  3.1-4,  and  the 
scattered  field  -  E^  is  shown  in  Figure  3.1-5.  These  results  are  for  the  value  X  =  D, 
the  plane  midway  through  the  slab  length.  Notice  that  for  parallel  polarization  a 
component  of  field  is  found  in  the  x  direction,  although  the  incident  field  has  no 
x  component. 

The  scattered  field  E^  computed  by  the  moment  method  has  an  exponential  decrease 
outside  the  slab  region.  For  a  TE  mode,  a  slab-guided  wave  decreases  as  exp(-v)  |y  ~|), 
where  v  =  0.704ir/inch.  For  a  TM  mode,  a  guided  wave  decreases  with  v  =  0.3367r/inch. 
The  exponential  decreases  are  graphed  in  Figure  3.1-6,  with  the  assumptions  that 
exponential  decrease  starts  at  y  =  0.1  inch  and  that  the  field  magnitudes  are  assumed 
correct  at  that  value  of  y. 

The  minima  in  |  ET|  occur  at  Y  =  0.6  in.  for  the  TE  case  and  Y  =  1.0  in.  for  the  TM  case 
(Figure  3.1-4).  These  locations  correspond  to  phase  values  of  tt  in  Es,  as  seen  in 
Figure  3.1-C,  and  the  zero  value  of  phase  of  E*. 

We  conclude  that  guided  waves  exist  and  produce  effects  at  distances  as  large  as  one 
wavelength  from  the  slab.  These  distances  are  an  appreciable  part  of  the  diameter  of 
small  missile  radomes. 

In  Figure  3.1-4,  the  phase  of  E^  varies  with  Y  more  rapidly  for  the  TE  wave  (or 
perpendicular  polarization)  than  for  the  TM  wave  (parallel  polarization).  This  result  is 
significant  for  radomes  and  possibly  for  guided  wave  integrated  antennas. 

3.2  TWO-LAYER  SLAB,  LENGTH  3.18  WAVELENGTHS 

Because  the  calculations  for  the  homogeneous  slab  showed  appreciable  guided  wave 
magnitude,  additional  calculations  were  done  for  a  thicker  slab.  This  slab  consisted  of 
two  layers:  one  layer  was  identical  to  that  of  the  homogeneous  slab  in  the  preceding 
section  (Table  3.1-1);  the  second  layer  had  the  same  thickness,  but  a  different  dielectric 
constant.  The  parameters  are  listed  in  Table  3.2-1.  The  slab  was  analyzed  to  evaluate 
how  a  second  layer  affects  guided  wave  excitation  and  decrease  with  distance  from  the 
slab.  The  wavelength  used  was  1.259  in. 

Table  3.2-1.  Parameters  of  Two- Layer  Slab 


Length: 

4  in. 

Thickness: 

layer  1 

0.125  in. 

layer  2 

0.125  in. 

Dielectric  Constant: 

layer  1 

2.6 

layer  2 

4.0 

Figure  3.2-1  shows  |e"T|  for  x  =  0,  the  plane  dividing  the  slab  symmetrically.  By 
comparing  Figures  3.2-1  and  3.1-4  for  parallel  polarization,  we  see  larger  maximum 
values  of  |E*|  for  the  thicker  slab.  In  contrast,  for  perpendicular  polarization  the 
maximum  value  of  |ET|  is  smaner  for  the  thicker  slab. 
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Figure  3.1-4.  Total  Field  on  Transverse  Plane  X  =  0,  for  4-Inch  Long  Slab: 
Thickness  0.125-in.;  Grazing  Incidence  Cell  Radii  0.0715-in.; 
Frequency  9.4  GHz;  Dielectric  Constant  2.6 
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Figure  3.1-5.  Scattered  Field  for  Figure  3.1-4 
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Figure  3.1-6.  Magnitude  of  Scattered  Field;  by  Moment  Method  (•); 

by  exp  (-v  Jy  -0.1|)- 

The  phase  of  ET  is  shown  in  Figure  3.2-2.  The  changes  of  phase  with  Y  have  opposite 
signs  for  the  two  polarizations.  If  the  slab  is  taken  to  be  a  model  of  one  side  of  an 
axially  symmetric  radome  with  grazing  incidence  on  that  side,  then  the  difference  in  the 
signs  of  the  spatial  derivative  of  phase  implies  opposite  signs  of  boresight  error. 

Figure  3.2-3  shows  |eS|,  and  Figure  3.2-4  shows  arg  E5*.  Significant  differences  exist 
for  the  two  polarizations. 

3.3  HOLLOW  WEDGE  WITH  WALLS  OF  LENGTH  3.18  WAVELENGTHS 

As  a  step  toward  analyzing  axially  symmetric  structures,  moment-method  calculations 
were  done  for  a  hollow  wedge  consisting  of  two  slabs  like  those  in  Section  3.2.  The 
parameters  of  the  wedge  are  listed  in  Table  3.3-1,  and  its  dimensions  are  shown  in  the 
coordinate  system  of  Figure  3.3-1.  The  wavelength  used  was  1.259  in. 

The  total  field  magnitude  |  j  inside  the  wedge  is  shown  in  Figure  3.3-2  for 
perpendicular  polarization  and  in  Figure  3.3-3  for  parallel  polarization.  The  figures  also 
show  |eT|  for  grazing  incidence  on  a  single  slab.  The  values  of  | E^l  for  grazing 
incidence  on  a  slab  are  very  close  to  those  for  the  wedge's  lower  wall,  which  is  also 
illuminated  at  grazing  incidence. 
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Figure  3.2-1.  Total  Field  for  Two- Layer  Slab,  Grazing  Incidence,  on  Plane  x  =  0 
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Figure  3.2-2.  Phase  of  E^  Corresponding  to  Figure  3.2-1 

For  axial  incidence,  Figure  3.3-4  shows  the  z  component  of  the  total  field  e£  in  the 
symmetric  transverse  plane  (x  =  2  in.)  for  perpendicular  polarization.  Figure  3.3-5  shows 
E^  for  parallel  polarization.  Notice  that  |E^|  is  not  zero;  this  is  becuase  the  wedge 
generates  a  field  component  in  the  direction  of  propagation  of  the  incident  wave. 

The  phase  curves  in  Figures  3.3-4  and  3.3-5  are  different  for  the  two  polarizations.  In 
addition,  the  phase  values  depend  on  the  direction  of  the  incident  wave.  The  dependence 
on  incidence  direction  results  in  boresight  error,  and  differences  in  the  phase  curves 
mean  that  boresight  error  depends  on  polarization. 
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Figure  3.2-4.  Phase  for  Figure  3.2-3 
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Table  3.3-1.  Parameters  of  Wedge 


Length  of  Sides 

4  in. 

Thickness  of  Walls 

0.125  in. 

Half  Angle 

16.5 

Dielectric  Constant: 

2.6 

To  quantify  these  effects  on  boresight  error  we  measured  the  phase-delay  difference 
between  two  points  symmetrically  located  about  the  wedge  axis  with  their  centers 
2.25  in.  (1.8  wavelength)  from  the  wedge  vertex.  Figure  3.3-6  shows  the  phase 
differences  between  the  pairs  of  points.  Notice  that  the  differences  for  perpendicular 
polarization  have  larger  magnitudes  than  those  for  parallel  polarization.  This  result  is 
significant  for  electronic  compensation  for  it  implies  the  need  for  polarization  sensing. 

In  the  calculations  for  parallel  polarization,  the  rectangular  components  Ex  and  Ey  were 
projected  onto  the  baseline  to  simulate  an  aperture  antenna. 


Figure  3.3-1.  Coordinates  for  Wedge 
3.4  AXIALLY  SYMMETRIC  SHELL 

The  shell's  outer  contour  is  shown  in  Figure  3.4-1.  It  is  a  curve  with  radius 
r  =  R(x  -L)1/2 

where  R  is  the  base  radius  and  x  is  the  axial  coordinate.  The  radome  was  fabricated 
from  a  resin,  with  dielectric  constant  2.1,  and  its  thickness  tapers  from  0.175  in.  at  the 
tip  to  0.130  in.  at  the  base. 
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Figure  3.3-2.  Total  Field;  Perpendicular  Polarization,  Inside  Wedge,  for 
Grazing  Incidence  on  Lower  Wall,  and  Total  Field  Inside  Slab. 
Number  of  cells:  wedge,  65;  slab,  33 


Figure  3.3-4.  Total  Field  at  x  =  2-in.  for  Wedge,  Axial  Incidence, 
Perpendicular  Polarization 

3.4.1  WAVEFRONT  MEASUREMENTS 

The  phase  and  intensity  irside  the  shell  were  measured  by  scanning  a  small  probe  antenna 
(aperture  l. 25-in.  square)  over  a  linear  path  with  approximately  plane  waves  incident 
externally.  The  receiver  was  a  network  analyzer,  and  a  waveguide  path  propagated  the 
coherent  reference  field.  The  radome  was  supported  by  an  absorber-covered  frame  that 
provided  rotation  about  a  vertical  axis  to  vary  the  incidence  direction  defined  in 
Figure  3.4-2  by  the  gimbal  angle  /3,  which  is  between  the  incident  wave  normal  and  the 
radome  axis.  The  radome  was  pivoted  about  a  point  on  the  radome  axis  1.5  in.  from  the 
base,  and  the  probe  aperture  scanned  a  horizontal  line  1.5  in.  closer  to  the  tip.  The 
frequency  used  was  10.000  GHz  (wavelength  =  1.180  in.). 

Measured  power  transmittance  and  phase  delay  are  shown  in  Figures  3.4-2  and  3.4-3  for 
horizontal  polarization;  this  scan  direction  and  polarization  are  the  E-plane. 
Figures  3.4-4  and  3.4-5  show  transmittance  and  phase  delay  for  vertical  polarization; 
this  case  is  referred  to  as  the  H-plane. 

Figure  3.4-2  shows  very  low  transmittance  values  for  horizontal  displacements  near  2  in. 
The  phase  delay  curves  in  Figure  3.4-3  have  overall  slopes  that  change  from  positive  to 
negative  as  /3  is  changed.  This  change  in  slope  implies  a  change  in  the  sign  of  boresight 
error  with  gimbal  angle  and  a  zero  value  at  approximately  6°.  In  contrast,  for  vertical 
polarization  the  slope  is  negative  for  most  gimbal  angles  although  for  larger  angles,  like 
16*,  the  phase  tilt  is  becoming  positive. 
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Figure  3.3-6.  Phase  Difference  Between  Two  Points; 

Separated  by  0.9  in.  at  x  =  2.25  in.  within  Wedge 

The  low  transmittance  values  near  the  wall  for  the  E-plane  are  surprising  because,  for 
flat  sheets,  parallel  polarization  has  high  transmittance.  The  low  values  suggest  guided 
waves,  as  do  the  values  that  exceed  0  dB. 

3.4.2  BORESIGHT  ERROR  MEASUREMENTS 

Boresight  error  was  estimated  by  measuring  the  phase  difference  between  a  pair  of  horn 
antennas.  The  antennas  were  connected  to  a  network  analyzer.  Each  horn  had  E-plane 
aperture  1.7  in.  and  H-plane  aperture  2.1  in.  The  antennas  were  spaced  by  2.1  in.  in  the 
H-plane.  Measurements  were  made  by  pivoting  the  radome  and  recording  phase 
difference.  The  antennas  were  used  in  two  angular  orientations,  rotated  about  the 
symmetry  axis  of  the  pair,  to  measure  E-  and  H-plane  errors;  incident  wave  polarization 
also  was  rotated. 

Figures  3.4-6  and  3.4-7  show  measured  phase  differences.  Note  the  sign  reversal  in  the 
E-plane  data  near  ±6*  gimbal  angle.  That  is,  the  phase  differences  for  13  less  than  6* 
have  signs  opposite  to  those  for  (3  exceeding  6*.  This  sign  reversal  does  not  occur  for  the 
H-plane  (Figure  3.4-7).  The  sign  reversal  in  the  E-plane  and  its  absence  in  the  H-plane 
are  reasonable  when  we  consider  the  probing  data  (see  the  discussion  in  subsection  3.4.1). 
Moreover  this  behavior  is  suggested  by  the  moment-method  results  in  Figure  3.3-6. 
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Figure  3.4-2.  Intensity  Transmittance  of  Axially  Symmetric  Shell. 
The  inset  shows  the  probing  path  as  a  dotted  line;  the  graphs 
are  labeled  with  values  of/3,  horizontal  polarization 
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3.4.3  WAVEFRONT  ANALYSIS 

This  section  describes  the  analysis  of  the  measured  wavefronts  into  two  separate  fields. 
The  procedure  starts  with  the  computation  of  the  wavefront  field  Ec  by  means  of  the 
surface  integration  method  (Figure  2-2).  The  second  step  is  to  subtract  Ec  from  the 
measured  values  Em.  The  difference  Em  -  Ec  is  then  analyzed  to  determine  its 
properties. 

The  details  of  computing  Ec  are  as  follows.  The  integration  is  confined  to  an  area 
defined  by  the  first  Fresnel  zone  which  is  defined  by  the  transverse  radius  r  =v/\l 
where  X  is  wavelength  and  i  is  the  distance  from  the  observation  point  (on  the  probing 
surface)  to  the  incident  wavefront  surface  containing  the  Fresnel  zone.  The  validity  of 
using  the  first  Fresnel  zone  has  been  discussed  in  another  report  (Reference  7).  The 
integration  mesh  size  is  not  very  critical  so  long  as  the  spacing  of  convergent  rays,  in  the 
Fresnel  zone,  is  approximately  a  half-wavelength.  The  value  of  i  is  a  critical  parameter 
because  it  influences  sampling,  therefore  Ec  was  calculated  for  i  equal  to  10,  20,  30,  40, 
and  50  wavelengths.  The  computed  values  of  Ec  were  compared  with  Em,  and  those  for 
I  equal  40X  agreed  most  closely  with  the  measured.  This  value  of  1  gives  an  illuminated 
area  on  the  radome  that  agrees  rather  well  with  that  found  earlier  for  a  shorter  radome 
that  had  the  same  base  diameter  as  the  shell,  but  was  only  16.5  in.  long. 

Figure  3.4-8  shows  samples  of  intensity  (|EC|2  or  |Em|2)  and  phase  (Ec  or  Em).  The 
discrepancies  in  phase  delay  are  largest  for  (the  transverse  coordinate  in  the  probing 
path),  approximately  equal  to  2  in.  The  region  in  which  the  probe  approaches  the  wall 
that  is  illuminated  at  approximately  grazing  incidence,  is  shown  in  the  inset.  The 
radome  offset  (gimbal)  angle  /3  was  15*.  The  discrepancies  of  phase  in  Figure  3.4-8  imply 
significant  discrepancies  in  boresight -error  calculations  that  omit  guided  waves. 

The  difference  between  Em  and  Ec  was  calculated  and  |Em  -  Ec|  and  arg  (Em  -  Ec) 
plotted  for /3  equal  15°  (Figure  3.4-9).  Also  plotted  are  two  curves  that  describe  the 
exponential  decrease  of  field  magnitude  for  a  slab-guided  wave.  Specifically,  the  graphs 
show  0.26  exp(-v|y-2|)  where  0.26  is  the  value  computed  for  |Em  -  Ec|  .  The  two 
exponential  curves  are  for  v  equal  0.75/inch  and  1.01/inch;  these  values  were  calculated 
for  infinite  slabs  with  thicknesses  0.125  in.  and  0.150,  respectively. 

The  agreement  between  |Em  -  Ec|  and  the  exponential  curves  supports  the  view  that 
the  radome  causes  guided  waves,  and  that  these  guided  waves  are  of  the  relatively 
simple  type  that  propagate  on  slabs.  The  magnitude  of  the  guided  wave  is  appreciable, 
but  of  more  importance  is  the  phase. 

The  phase  variation  of  Em  -  Ec  (Figure  3.4-9)  can  also  be  explained  in  terms  of  slab- 
guided  waves.  If  we  assume  a  slab-guided  wave,  the  equiphase  surface  in  Figure  3.4-10 
describes  the  phase  surface.  This  phase  surface  crosses  the  probing  path  at  angle  of 
approximately  10*,  not  15*,  because  the  radome  has  axial  curvature.  The  phase  variation 
caused  by  the  assumed  guided  wave  is  then  (kg  Ay  a  sin  6  ),  where  6  is  the  angle 
between  the  wavefront  and  the  probing  path,  is  displacement,  and  kg  is  the 

propagation  constant  of  the  slab-guided  wave.  For  thickness  =  0.125  in.,  kg  is 
1.720n-/inch,  and  for  thickness  =  0.150  in.,  kg  is  1.75lir/inch.  A  graph  of  phase  change 
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Figure  3.4-8.  Intensity  and  Phase  Axially  Symmetric  Shell  for  /?  =  15*;  E-Plane. 
Measured  with  1.25-in.  square  antenna  (x).  Calculated  with  surface 
integration  (•)  for  I  =  40A;  frequency  10  GHz. 
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Figure  3.4-9.  The  Magnitude  and  Phase  of  Em  -  Ec  (•).  Computed 

values  of  0.260  exp  (-v  |  yA  -2  | );  for  thickness  0.125  in.  ( - ); 

for  thickness  0.150  in.  ( - ). 


(kg  y a  sin  9  )  is  superimposed  on  the  phase  values  of  Em  -  Ec  in  Figure  3.4-9;  the 
location  of  the  line  on  the  ordinate  scale  is  arbitrary:  the  slope  is  significant. 

The  comparison  of  computed  and  measured  phase  in  Figure  3.4-9  further  confirms  the 
view  that  the  field  near  the  wall  of  the  shell  contains  a  constituent  that  is  a  guided 
wave,  and  the  constituent  is  a  slab-guided  wave. 

This  model  omits  the  opposite  radome  wall  and  a  backward  propagating  guided  wave  as 
an  approximation.  Furthermore,  the  model  does  not  acknowledge  the  dependence  of 
guided  wave  excitation  on  gimbal  angle;  this  dependence  for  slabs  has  been  established  in 
our  earlier  work  with  the  moment  method.  Although  Em  -  Ec  has  been  calculated  for 
additional  gimbal  angles,  we  omit  the  results  for  brevity. 

It  is  useful  to  state  that  v  is  larger  for  perpendicular  polarization  than  for  parallel. 
Consequently,  the  exponential  function  decreases  less  rapidly  for  parallel  polarization. 
We  expect  guided  wave  effects  to  be  more  pronounced  in  E-plane  data. 

In  addition,  the  gimbal  angle  /3  plays  a  role.  For  fi  as  in  Figure  3.4-10,  the  phase  of 
Em  -  Ec  decreases  with  YA  (Figure  3.4-9)  because  the  points  on  the  probing  path  move 
nearer  the  source  of  the  wave  as  Y^  increases.  For  smaller  offset  angles,  the  phase 
surface  tilts  so  that  /3  is  negative  and  phase  increases  Y^.  This  behavior  explains  the 
sign  reversal  in  the  boresight  error  curve  of  Figure  3.4-6,  which  is  for  the  E-plane.  For 
the  H-plane  the  value  of  v  is  2.21/inch  for  thickness  0.125  in.;  the  guided  wave  is 
strongly  attenuated. 


Figure  3.4-10.  A  Suggestion  of  a  Guided  Wave  on  a  Slab  that  Approximates 
the  Side  of  the  Radome  Near  Y^  and  Equals  2-ln.  in  Figure  3.4-6 
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4.  HOLLOW  FINITE  CYLINDER 


Fields  near  a  finite  dielectric  cylinder  were  computed  using  the  new  version  of  a  moment 
method  theory  outlined  in  Section  2.  This  new  version  includes  an  axial  component,  as 
well  as  transverse  components,  of  field  near  the  scatterer.  The  newer  theory  gives 
better  agreement  with  measured  nearfields  than  did  an  older  theory  that  omitted  the 
axial  component.  Polarization  dependence  is  correctly  computed. 

4.1  COMPUTATIONS 

The  fields  near  a  hollow  finite  cylinder  were  computed  for  axial  incidence.  The 
dielectric  constant  was  2.6,  outside  diameter  was  1.25  in.,  thickness  was  0.065  in.,  and 
length  was  0.42  in. 

To  test  the  effect  of  cell  size  upon  the  fields,  two  cylinders  like  those  in  Figure  4.1-1 
were  used  to  evaluate  the  equations  of  the  preceding  section.  Each  cylinder  was  divided 
into  two  rings.  The  first  cylinder  had  16  cells  per  ring;  N  was  32  for  the  entire  cylinder. 
The  second  cylinder  had  32  cells  per  ring;  N  was  64.  The  external  field  at  1.26  in.  from 
the  center  of  the  ring  was  computed  for  cylinder  one  (Figure  4.1-2). 


Figure  4.1-1.  Subdivision  of  Hollow  Cylinder  into  Rings.  An  observation  point 
has  coordinates  <x',  y',  z');  it  may  be  inside  or  outside  the  cylinder. 

Calculations  for  the  second  cylinder  (N  =  64)  are  shown  in  Figures  4.1-3  and  4.1-4  for 
distance  1.26  in.,  and  in  Figures  4.1-5  and  4.1-6  for  distance  0.63  in.  The  calculated 
values  are  for  the  equations  in  Section  2,  with  Ez  #  0,  and  for  the  equations  of 
Reference  3,  with  Ez  =  0.  The  inclusion  of  Ez  and  the  value  of  N  influence  the  results. 

4.2  DISCUSSION 

Figures  4.1-2  through  4.1-6  show  measured  as  well  as  computed  values.  The  measure¬ 
ments  were  made  by  scanning  a  half-wave  dipole  probe  on  horizontal  and  vertical  lines 
(H-  and  E-planes). 
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Figure  4.1-2.  Total  Field  at  z  =  1.26  in.  for  0.42  in.  Long  Ring.  Measured  (x) 
calculated,  N  =  64,  Ez  #  0(  •  >;  calculated,  N  =  32,  Ez  **  0(+) 
Wavelength:  1.26  in.  dielectric  constant  2.6;  diameter  1.25  in. 


Figure  4.1-3.  Total  Field  at  z  =  1.26  in.  for  0.42  in.  Long  Ring.  Measured  (x); 

calculated,  N  =  64,  Ez  =  0  (0);  calculated  N  =  64,  Ez#0  (•).  Wave¬ 
length:  1.26  in.  Dielectric  constant:  2.6  diameter  1.25  in.  Calcula¬ 
tions:  32  cells/axial  subdivision 
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Figure  4.1-4.  Total  Field  at  z  -  1.26  in.  for  0.42  in.  Long  Ring.  Measured  (x); 

calculated,  N  =  64,  Ez  =  0  (0);  calculated  N  =  64,  Ez#0(»). 
Wavelength:  1.26  in.  Dielectric  constant  2.6,  diameter  1.25  in, 
Calculations:  32  cells/axial/subdivision 
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Figure  4.1-5.  Total  Field  at  z  =  0.63  in.  for  0.42  in.  Long  Ring.  Measured  (x); 

calculated,  N  =  64,  Ez  =  0  (0);  calculated  N  =  64,  Ez*0  (•). 
Wavelength:  1.26  in.  Dielectric  constant  2.6,  diameter  1.25  in. 
Calculations:  32  cells/axial/subdivision 
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Figure  4.1-6.  Total  Field  at  z  =  0.63  in.  for  0.42  in.  Long  Ring.  Measured  (x); 

calculated,  N  =  64,  Ez  =  0  (0);  calculated  N  =  64,  Ez#0  (•>. 

Wavelength:  1.26  in.  Dielectric  constant  2.6,  diameter  1.25  in. 
Calculations:  32  cells/axial/subdivision 


Figure  4.1-2  shows  that  the  discrepancies  between  measured  and  computed  values  are 
smaller  for  N  =  64  than  for  N  =  32. 

Figures  4.1-3  and  4.1-4  show  discrepancies  are  smaller  for  the  equations  with  Ez?*0  than 
they  are  for  Ez  =  0.  The  same  conclusion  follows  from  Figures  4.1-5  and  4.1-6. 

Accuracy  is  good  in  Figures  4.1-3  and  4.1-4  where  the  distance  between  the  probe  and 
scatterer  is  a  wavelength.  Notice  that  the  spatial  dependence  in  the  E-plane 
(Figure  4.1-3)  differs  from  that  in  the  H-plane  (Figure  4.1-4).  Accuracy  is  not  as  good 
when  the  probe-scat  ter  er  distance  is  half  a  wavelength  (Figures  4.1-5  and  4.1-6). 
Reference  5  showed  that  on-axis  discrepancies  were  reduced  by  introducing  a  probe 
correction. 


5.  HOLLOW  CONES 


This  section  presents  data  for  measured  and  computed  fields  for  a  hollow  cone  for  the 
two  versions  of  the  moment  method  theory  (Ez  =  0,  and  Ez  #  0)  outlined  in  Section  2. 
The  results  are  for  axial  and  nonaxial  incidence. 

It  was  shown  in  the  preceding  sections  that  the  shapes  and  sizes  of  the  cells  into  which  a 
dielectric  structure  is  decomposed  influence  computations.  Accuracy  depends  on  cell 
size  because  of  the  assumption  that  the  field  in  a  cell  is  constant;  therefore,  smaller 
cells  are  desirable.  However,  computation  increases  with  decreasing  cell  size,  and 
eventually  cell  sizes  can  be  made  so  small  that  the  storage  capacity  of  computers  may 
be  inadequate  to  solve  the  large  set  of  simultaneous  equations.  In  addition,  as  the  cell 
size  is  decreased,  the  computed  fields  might  not  converge  to  a  reasonably  fixed 
distribution  of  values.  Finally,  cell  shape  should  be  chosen  to  exploit  those  directions  in 
which  the  field  varies  slowly;  of  course,  the  cells  should  fit  the  scatterer's  geometry. 

With  these  limitations  in  mind,  we  have  chosen  for  hollow  cones  a  division  that  is 
described  in  Figures  2.2-1,  2.2-2,  5-1,  and  5-2.  Figure  5-1  is  a  picture  of  an 
experimental  model  that  consists  of  15  rings  of  Plexiglass;  another  version  with  13  rings 
was  also  utilized.  Figure  5-2  gives  dimensions  for  the  13-ring  cone. 


Figure  5-1.  Hollow  Cone  Composed  of  Rings.  Half  angle  is  14*. 
Dielectric  constant  was  2.6. 
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Figure  5-2.  Dimensions  of  Rings.  Cell  lengths  2  8z  were  0.215  in.  Thickness  25p  were 
0.065  in.  This  figure  shows  only  five  of  the  rings,  but  the  models  had  either 
13  or  15.  Overall  dimensions  in  this  figure  are  for  a  cone  with  13  rings. 

5.1  COMPUTATIONS  AND  MEASUREMENTS  FOR  HOLLOW  CONE 

This  section  describes  calculations  and  measurements  for  a  13-ring  cone  with  the 
dimensions  in  Figure  5-2.  Results  are  given  for  theory  with  Ez  #  0,  and  Ez  =  0,  for  axial 
and  nonaxial  incidence. 

Phase  and  intensity  were  measured  with  a  half-wave  dipole  antenna  in  the  region 
bounded  by  the  cone  and  over  planes  behind  the  cone.  The  receiver  was  a  network 
analyzer;  its  reference  field  propagated  through  a  waveguide.  The  transmitting  and 
receiving  antennas  were  spaced  by  100  inches.  The  cone  was  suspended  by  nylon  threads 
in  an  anechoic  chamber. 

The  procedure  for  positioning  the  probe  relative  to  the  cone  was  to  start  with  the  probe 
at  a  specific  axial  distance  from  the  base  and  to  scan  the  probe  laterally.  Axial  position 
was  changed  and  scans  repeated.  The  cone  was  stationary.  To  obtain  data  for  the 
incident  field  the  cone  was  removed.  This  procedure  was  more  accurate  than  an  earlier 
procedure  in  which  the  cone  was  removed  or  replaced  between  scans  and  axial  probe 
displacements. 

Figure  5.1-1  shows  computed  values  of  |E|2,  the  intensity  corresponding  to  the  total 
internal  field,  for  axial  incidence,  on  two  cone  generators.  One  generator  is  in  the  x-z 
plane,  called  the  E-plane  because  it  contains  the  electric  field.  The  other  generator  is  in 
the  y-z  plane,  called  the  H-plan?.  Note,  that  the  field  in  the  H-plane  is  larger  than  that 
in  the  E-plane.  The  results  differ  somewhat  for  the  two  13-ring  cases.  In  one  case  each 
ring  contained  12  cells.  In  the  other  case,  the  six  smallest  rings  had  eight  cells  each,  and 
the  seven  largest  rings  had  16  cells  each.  The  main  differences  are  near  the  tip,  within 
the  smaller  cells.  The  smallest  rings  have  thicknesses  approximately  equal  to  their  radii, 
a  condition  which  violates  an  assumption  in  the  derivation  of  the  working  formulas. 
Therefore,  calculations  were  made  by  representing  each  of  the  three  smallest  rings  of 
the  13-ring  cone  by  a  pair  of  concentric  rings  and  retaining  the  10  remaining 
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Figure  5.1-1.  Computed  |Ej  2  Inside  Conical  Shell.  The  inset  shows  the  loci  of  points. 

E  refers  to  points  in  the  x-z  plane,  the  E-plane;  H  refers  to  y-z  plane, 
the  H-plane.  Computed  for  Ez  =  0. 

rings  as  before.  Figure  5.1-2  shows  the  results.  It  is  clear  that  the  field  depends  on 
axial  position. 

To  evaluate  the  significance  of  the  difference  between  the  internal  fields  for  the  two 
subdivisions,  we  compared  the  respective  external  fields  computed  for  the  two  sets  of 
internal  fields.  Figure  5.1-3  shows  the  measured  and  computed  fields  for  both  sets  at 
one-half  wavelength  behind  the  largest  ring.  Near  the  main  maxima  for  phase  and 


intensity,  the  results  agree  rather  well,  but  considerable  differences  exist  in  intensity  at 
the  deepest  minimum  and  the  secondary  maxima.  At  the  intensity  maximum,  the 
discrepancy  between  the  measured  value  and  that  computed  with  13  rings  is  approxi¬ 
mately  8%;  the  discrepancy  is  approximately  15%  for  16  rings.  The  intensity  dis¬ 
crepancies  at  the  first  minimum  and  subsidiary  maximum  favor  neither  subdivision. 

For  axial  incidence.  Figure  5.1-4  shows  the  measured  and  computed  values  of  phase 
delay  and  |EX|  2,  the  intensity  corresponding  to  the  x  component  of  field.  The  field  was 
evaluated  at  z  =  0.86  in.,  which  is  the  plane  of  the  center  of  the  12th  ring,  that  next  to 
the  largest.  The  results  are  for  both  the  E  and  H  planes,  and  for  Ez  =  0  and  Ez  *0. 

For  nonaxial  incidence  (z  =  0.86  in.)  Figure  5.1-5  shows  A$and  |EX|^  for  the  incident 
wave  normal  at  14.9*  to  the  cone  axis.  Incidence  is  nearly  grazing  on  one  side  of  the 
cone.  For  axial  incidence,  Figures  5.1-6  and  5.1-7  show  the  computed  results  for  Ez  #  0 
at  z  =  1.80  in.  (half  wavelength  behind  the  largest  ring)  for  the  E  and  H  planes. 

For  the  plane  one  half  wavelength  behind  the  largest  ring,  Figures  5.1-8  and  5.1-9  show 
phase  and  intensity  (normalized  to  the  incident  field)  for  the  E-  and  H-planes  for  axial 
incidence.  For  off-axis  incidence  (14.9*  between  the  cone  axis  and  incident  wave 
normal).  Figures  5.1-10  and  5.1-11  show  phase  and  intensity  for  the  E-  and  H-planes, 
respectively,  for  off-axis  incidence. 

Figure  5.1-12  gives  measured  values  with  the  probe  antenna  on  the  x-axis  at  distance 
half  wavelength  behind  the  cone,  and  Figure  5.1-13  gives  measured  values  at  the  same 
distance  for  the  probe  on  the  y-axis.  Two  sets  of  values  are  shown:  one  set  is 
measurements  taken  in  August  1982;  the  other  set,  in  August  1980.  The  intensity  data 
are  quite  similar,  but  the  phase  curves  differ.  This  difference  is  approximately  constant. 
'Me  attribute  it  to  imprecision  in  positioning  the  cone.  The  quantity  a<£  is  determined  by 
comparing  phases  measured  with  and  without  the  cone.  In  the  1980  set  of  measurements, 
the  probe  scanned  lines  at  specific  axial  positions  with  the  cone  present;  then  the  cone 
was  removed  and  the  probe  repositioned.  In  the  1982  set,  the  probe  was  not  moved 
axially  between  measurements  with  and  without  the  cone. 

5.2  DISCUSSION 

Preceding  sections  showed  that  accuracy  was  greater  for  calculations  that  included  both 
the  axial  and  transverse  components  than  for  those  with  only  the  transverse  component. 

There  was  greater  accuracy  for  the  probe  positioned  a  half  wavelength  behind  the  base 
of  the  cone,  but  for  measurements  inside  the  region  bounded  by  the  cone,  in  the  plane  of 
the  12th  ring,  the  relative  accuracy  was  less  obvious.  The  effects  of  the  probe  may  be 
greater  when  it  is  nearer  the  cone.  Figure  5.1-4  shows  a  strong  polarization 
dependence:  |  E  |  and  A4  vary  more  for  the  H-plane  than  for  the  E-plane.  The  values 
of  |  E  |  are  larger  for  the  H-plane,  as  expected  from  Figure  5.1-1. 

Figure  5.1-4  also  shows  that  discrepancies  in  A are  smaller  for  the  theory  that  has 
Ez#0  than  they  are  for  the  theory  with  Ez  =  0.  For  Ez*0,  phase  discrepancies  range 
from  0  to  15%;  for  Ez  =  0,  discrepancies  range  from  2  to  25%.  However,  intensity 
discrepancies  are  greater  for  Ez#0,  ranging  from  3  to  20%,  and  they  equal  10%  for 
F.z  =  0. 
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Figure  5.1-4.  Phase  (A*)  and  Intensity  (|E^|  2)  Measured  (X)  and  Computed  in  Plane 
of  12th  Ring  (at  z  =  0.86  in.)  for  13-Ring  Cone.  The  electric  field  was 
horizontal,  in  the  x-z  plane.  The  H-plane  is  the  y-z  plane  (x  =  0)  and  the 
E-plane  is  the  plane  y  =  0,  as  shown  in  the  sketches  above.  The  probe 
was  a  half-wave  dipole.  Calculations  were  done  with  eight  cells  per  ring. 
For  Ez  =  0  (o);  for  Ez#0  (•). 
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Figure  5.1-5.  Phase  (A 4)  and  Intensity  (|E|2)  Measured  and  Computed  in  Plane  of 
12th  Ring  (at  z  -  0.86  in.)  of  13-Ring  Cone,  for  14.9*  Off-axis  Angle  of 
Incident  Wave.  The  E-field  was  horizontal;  E-  and  H-planes  are  defined 
in  Figure  5.1-4.  The  probe  was  a  half-wave  dipole.  Computed  with  eight 
cells  per  ring.  For  Ez  =  0  (o);  for  Ez#0(>). 


Figure  5.1-6.  Phase  (A*)  and  Intensity  (|Ex|2)  at  z  =  1.81  in.,  a  Half  Wavelength 
Behind  Largest  Ring  of  13-Ring  Cone.  Computed  (•),  with  EZf*0,  with 
eight  cells  per  ring;  measured  (x);  wavelength  1.26  in.  Incidence  was 
axial.  Data  are  for  the  E-plane. 
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H-Plane  Phase  (A#)  and  Intensity  (|EX|2)  z  =  1.81  in.,  a  Half 
Wavelength  Behind  Largest  Ring  of  13-Ring  Cone.  Computed  (•),  with 
Ez9*0,  with  eight  cells  per  ring;  measured  (x);  wavelength  1.26  in. 
Incidence  was  axial. 
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>ase  and  Intensity  Computed  for  z  =  1.81  in.,  One  Half  Wavelength 
(hind  Largest  Ring.  Computations  were  done  for  13  rings,  with  eight 
Us/ring.  For  Ez  =  0(0;  for  Ez#0  (+).  E- Plane. 
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e  and  Intensity  at  z  =  1.81  in.  for  13-Ring  Cone.  Calculated 
d  (x).  The  incident  wave  normal  was  14.9’  off  axis. 
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Figure  5.1-12.  Measured  A$and  I E  | 2 
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Figure  5.1-13.  As  in  Figure  5.1-12  but  for  the  H-Plane 
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For  nonaxial  incidence.  Figure  5.1-5  shows  A$  and  |EX|^  in  the  plane  of  the  12th  ring. 
Discrepancies  are  small,  and  the  curves  have  correct  shapes.  For  the  E-plane,  phase 
discrepancies  are  between  2  and  8%.  For  the  H-plane,  the  discrepancy  for  the  axial 
position  is  30%,  and  it  is  20%  for  off-axis  points;  these  discrepancies  equal  approxi¬ 
mately  6'  in  phase.  In  Figure  5.1-4  the  measured  and  computed  graphs  of  A<&  for  the 
E-plane  have  opposite  shapes,  but  discrepancies  are  still  small. 

Figures  5.1-8  and  5.1-9  show  A*  and  |EX|^  one  half  wavelength  behind  the  cone,  for 
axial  incidence. 

The  graphs  for  the  E-plane  and  the  H-plane  have  different  shapes.  The  secondary 
maxima  for  the  H-plane  are  higher  than  those  for  the  E-plane.  Discrepancies  of  A$  for 
x  equal  zero  are  approximately  8*,  or  25%.  Discrepancies  in  phase  are  smaller  at  off- 
axis  positions,  but  the  percentage  discrepancies  (ratio  of  discrepancy  to  measured  value) 
are  large  because  the  values  are  small.  The  shapes  of  the  theoretical  and  measured 
curves  agree  well  as  do  locations  of  extrema. 

Figures  5.1-10  and  5.1-11  show  A*  and  Iex|  ^  one  half  wavelength  behind  the  cone,  for 
off-axis  incidence.  Again  the  E-plane  and  H-plane  curves  have  different  shapes.  The 
H-plane  intensity  maxima  at  off-axis  positions  are  as  high  as  that  for  the  axial  position 
(y  equals  zero);  however,  the  E-plane  intensity  maximum  at  x  equal  0.5  in.  is  much 
higher  than  the  subordinary  maxima.  The  phase  and  intensity  plots  are  symmetric  in  the 
H-plane,  but  in  the  E-plane  the  off-axis  incidence  skews  the  data.  The  maximum  is 
shifted  from  the  axis  toward  the  direction  of  the  advancing  wave  normal.  The  theory 
agrees  well  with  the  measurements  although  at  the  discrepancies  phase  maximum  in 
Figure  5.1-10  are  30%.  In  Figure  5.1-11  the  phase  discrepancy  at  the  maximum  is  17%. 
These  large  values  may  be  caused  in  part  by  the  probe  averaging  the  field. 

Comparison  of  the  phase  data  in  Figures  5.1-8  and  5.1-12  shows  smaller  discrepancies 
for  Ez#0  than  for  Ez*0.  The  intensity  discrepancies  are  smaller  for  Ez#0  near  the 
maxima  and  minima,  but  the  computed  curves  intersect  near  inflection  points. 

Comparison  of  Figures  5.1-9  and  5.1-13  shows  smaller  phase  and  intensity  discrepancies 
for  Ezf*0  than  for  Ez  =  0. 

The  theories  give  the  polarization  difference,  which  appears  in  the  differences  between 
the  shapes  of  the  phase  curves  in  Figures  5.1-8  and  5.1-12  on  the  one  hand  and  those  in 
Figures  5.1-9  and  5.1-13  on  the  other. 

Discrepancies  seem  reasonable  in  light  of  the  fact  that  the  calculations  used  only  eight 
cells  per  ring.  The  inclusion  of  Ez  seems  to  improve  the  theory. 
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6. 


CONCLUSIONS 


This  project  produced  several  conclusions  about  electromagnetic  waves  near  dielectric 
slabs,  shells,  cylinders,  and  shells.  These  structures  are  models  of  radomes,  which 
produce  boresight  error  that  degrades  the  operation  of  avionic  systems  and  guidance 
systems. 

The  report  gives  specific  experimental  and  numerical  evidence  for  guided  waves  on  slabs 
and  axially  symmetric  shells.  Such  waves  are  omitted  in  conventional  radome  analysis. 
However,  they  have  significant  influence  on  broadband  systems.  They  influence 
transmittance  and  boresight  error,  and  are  polarization  dependent. 

The  report  describes  an  extension  of  the  moment  method.  Although  our  theory  is  based 
on  Richmond’s  vector  and  scalar  potential  approach,  we  did  analysis  for  a  new 
decomposition  of  curved  shells  and  derived  closed  form  expressions  for  matrix  elements 
without  appealing  to  regular  shapes  or  infinite  sizes  of  cells.  Moreover,  the  theory 
treated  polarization  by  including  both  axial  and  radial  field  components.  The  addition  of 
the  axial  component  complicated  the  theory. 

Experiments  with  cylinders  and  cones  tested  the  accuracy  of  nearfield  calculations. 
Discrepancies  between  experiment  and  calculation  were  the  criteria.  Small 
discrepancies  were  produced  by  the  theory  that  included  both  axial  and  radial 
components  even  for  off-axis  incidence.  Discrepancies  were  appreciable  for  a  version  of 
the  theory  that  had  only  two  radial  fields. 
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